Liquid-Gas phase transition in Bose-Einstein Condensates with time evolution 
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We study the effects of a repulsive three-body interaction on a system of trapped ultra-cold atoms 
in Bose-Einstein condensed state. The stationary solutions of the corresponding s— wave non-linear 
Schrodinger equation suggest a scenario of first-order liquid-gas phase transition in the condensed 
state up to a critical strength of the effective three-body force. The time evolution of the condensate 
with feeding process and three-body recombination losses has a new characteristic pattern. Also, 
the decay time of the dense (liquid) phase is longer than expected due to strong oscillations of the 
mean-square-radius. 
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The experimental evidences of Bose-Einstein conden- 
sation (BEG) in magnetically trapped weakly interacting 
atoms |^-|^ brought a considerable support to the theo- 
retical research on bosonic condensation. The nature of 
the effective atom-atom interaction determines the stabil- 
ity of the condensed state: the two-body pseudopoten- 
tial is repulsive for a positive s— wave atom-atom scat- 
tering length and it is attractive for a negative scattering 
length [Q. The ultra-cold trapped atoms with repulsive 
two-body interaction undergoes a Bose-Einstein phase- 
transition to a stable condensed state, in a number of 
cases found experimentally, as for ^^Rb for ^'^Na |^] 
and ^Li [||. However, a condensed state of atoms with 
negative s— wave atom-atom scattering length would be 
unstable for a large number of atoms |^,^ . It was indeed 
observed in the ^Li gas for which the s— wave scat- 
tering length is a = (—14.5 ± 0.4) A, that the number 
of allowed atoms in the condensed state was limited to a 
maximum value between 650 and 1300, which is consis- 
tent with the mean- field prediction 

From a theoretical approach, the addition of a repul- 
sive three-body interaction can extend considerably the 
region of stability for a condensate even for a very weak 
three-body force |l^. As one can observe from Refs. 
both signs for the three-body interaction are, in prin- 
ciple, allowed. However, in the present study we only 
consider the case of a repulsive three-body elastic inter- 
action together with an attractive two body interaction. 
We will show that, due to the repulsive three-body force, 
new physical aspects appears in the time evolution of the 
condensate. In respect to the static situation, it was sug- 
gested that, for a large number of bosons the three-body 
repulsion can overcome the two-body attraction, and a 
stable condensate will appear in the trap Singh 
and Rokhsar |l^ have also observed that above a critical 
value the only local minimum is a dense gas state, where 
the neglect of three-body collisions fails. 

In this work, using the mean-field approximation, we 
develop the scenario of collapse, which includes two as- 
pects of three-body interaction, that is recombination 
and replusive mean-field interaction. We begin by in- 
vestigating the competition between the leading term of 



an attractive two-body interaction, which is originated 
from a negative two-atom s— wave scattering length, and 
a repulsive three-body interaction, which can happen in 
the Efimov limit |l^, when \a\ oo. (The physics of 
three- atoms in the Efimov limit is discussed in Refs. p2|). 
We first consider the stationary solutions of the cor- 
responding extension of the Ginzburg-Pitaevskii-Gross 
(GPG) nonlinear Schrodinger equation (NLSE), for 
fixed number of particles, without dissipative terms, ex- 
tending an analysis previously reported in Refs. |l^. The 
liquid-gas phase transition in the condensate, suggested 
in , was confirmed by a more detailed analysis in the 
present stationary calculations. Then, the time evolu- 
tion of the feeding process of the condensate by an ex- 
ternal source is obtained by solving the time-dependent 
NLSE with repulsive three-body interaction (given by 
(73 > 0) and dissipation due to three-body recombina- 
tion processes. The dramatic collapse and the conse- 
quent atom loss that happens at the critical number of 
atoms (when — 0) [l^l is softened by the addition of 
the three-body repulsive force. The decay time of the 
liquid-phase is also unexpectedly long, when compared 
with the decay time that occurs for = 0, which gives a 
clue about the possible observation of three-body inter- 
action effects. Our results pointed out that the mean- 
square-radius is an important observable to be analyzed 
experimentally to study the dynamics of the growth and 
collapse of the condensate In the present study, in 

order to enphasize the real part of the three-body inter- 
action, we choose g^ significantly larger than the magni- 
tude of the dissipative term; although, in general, they 
are expected to be of the same order. 

The NLSE, which describes the condensed wave- 
function in the mean-field approximation, after consid- 
ering the two-body attractive and three-body repulsive 
effective interaction, is variationally obtained from the 
corresponding effective Lagrangian (see Gammal et al. 
in p3[). By considering a stationary solution, '^{r,t) = 
^-ifiTjh ^(jf^ where ji is the chemical potential and ipir) 
is normalized to the number of atoms N . By rescaling 
the NLSE for the s— wave solution, we obtain 
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for a < 0, where a; = yjimuj/h r and </'(x) = 
•\/87r|a| r'4!{r). The dimensionless parameters, related 
to the chemical potential and the three-body strength 
are, respectively, given by /3 = ji/fiuj and = 
X^hum? / (AttH'^ a)"^ . The normalization for 4>{x) reads 
/p dx\(j){x)\'^ ~ n where the reduced number n is re- 
lated to by n = 2N\a\y/2mLL>/h. The boundary con- 
ditions Q in Eq.(|l|) are given by 0(0) = and 
-> Cexp(-a;V4 + [/? - 5] ln(a;)) when x ^ 00. @ 
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FIG. 1. Central density pc, total energy E, chemical po- 
tential fi, and average square radius (f"^), as functions of 
the reduced number of atoms n. The three-body strength 
(?3 are given in the upper frame. The corresponding units 
are: (mw//i)/(47r|a|) for pc, {NTiu))/n for E, huj for /i, and 
n/[2mu)) for (r^). 

In Fig. 1, considering several values of 53 (0, 0.012, 
0.015, 0.0183 and 0.02), using exact numerical calcula- 
tions, we present the evolution of some relevant physical 
quantities -E, pc and (r^) as functions of the reduced 
number of atoms n. For — Q, our calculation repro- 
duces the result presented in Ref. |||,|l^, with the max- 
imum number of atoms limited by Umax ~ 1.62 (n is 
equal to |C3f I of Ref. §). 

As shown in the figure, for < 33 < 0.0183, the 
density pc, the chemical potential /i and the root-mean- 
squared radius (r^) present back bendings typical of a 



first order phase transition. For each g^, the transition 
point given by the crossing point in the E versus n cor- 
responds to a Maxwell construction in the diagram of fi 
versus n. At this point an equilibrated condensate should 
undergo a phase transition from the branch extending to 
small n to the branch extending to large n. The system 
should never explore the back bending part of the dia- 
gram because as we have seen in Fig. 1 it is an unsta- 
ble extremum of the energy. From this figure it is clear 
that the first branch is associated with large radii, small 
densities and positive chemical potentials while the sec- 
ond branch presents a more compact configuration with 
a smaller radius a larger density and a negative chemi- 
cal potential. This justify the term gas for the first one 
and liquid for the second one. However we want to stress 
that both solutions are quantum fluids. With 33 = 0.012 
the gas phase happens for n < 1.64 and the liquid phase 
for n > 1.64. For 53 > 0.0183 all the presented curves 
are well behaved and a single fluid phase is observed. At 
ga « 0.0183 and n « 1.8, the stable, metastable and un- 
stable solutions come to be the same. This corresponds 
to a critical point associated with a second order phase 
transition. At this point the derivatives of p,, pc and (r^^ 
as a function of n all diverge. We also checked that cal- 
culations with the variational expression of (r^), pc and 
fi are in good agreement with the ones depicted in Fig. 
1. 

In the lower frame of Fig. 2, wc show the phase bound- 
ary separating the two phases in the plane defined by n 
and (73 and the critical point at n « 1.8 and 53 « 0.0183. 
In the upper frame, we show the boundary of the forbid- 
den region in the central density versus 33 diagram. 
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FIG. 2. Phase diagram of the Bose condensate. Central 
density pc in units of (mu) /h)/{4:TT\a\). 

The main physical characteristic of the repulsive three- 
body force is to prevent the collapse of the condensate 
for the particle number above the critical number found 
with only two-body attractive interaction. The three- 
body repulsive potential tends to overcome the attrac- 
tion of the two-body potential at short distances, as de- 
scribed by Eq. (3), as the repulsive interaction grows 
inversely to x^, while the two-body potential is propor- 
tional to x^^. Thus, the implosive force that shrinks the 
condensate at the critical number is compensated by the 
repulsive three-body force. The time evolution of the 
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growth and collapse of the condensate with attractive in- 
teractions should be qualitatively modified by the 
presence of the repulsive three-body force. The three- 
body recombination effect jl^, which "burns" partially 
the condensed state should be taken into account to de- 
scribe quantitatively the dynamics of the condensate. In 
the case of only two-body attractive interaction, as ob- 
served by Kagan et al. [Q, by considering the feeding of 
the condensate from the nonequilibrium thermal cloud, 
the time evolution is dominated by a sequence of growth 
and collapse of the trapped condensate. The collapse 
occurs when the number of atoms in the condensate ex- 
ceeds the critical number N^, and it is followed by an 
expansion after the atoms in the high density region of 
the wave-function are lost due to three-body recombina- 
tion processes and consequently the average attractive 
potential from the two-body force is weakened. It is also 
noticed in Ref. that the condensate time evolution 
is dominated by an oscillatory mode of frequency ^ w; 
and, as time grows and N reaches a value > Nc, a huge 
compression takes place to implode the condensate. The 
repulsion given by the three-body force will dynamically 
affect the compression of the condensate weakening the 
implosive force and allowing more atoms to survive at 
high densities. 

In order to quantitatively study the above features 
with repulsive three-body interaction, we consider the 
time-dependent non-linear Schrodinger equation corre- 
sponding to Eq. including three-body recombination 
effects (with an intensity parameter 2^) and an imaginary 
linear term corresponding to the feeding of the conden- 
sate (with intensity parameter 7): 



d2 x2 |$|2 ^ l$|4 i^' 



where $ = ^(x, r) and r = cut. For the parameters ^ and 
7 we are using the same notation as given in Ref. | p^ . 

In Fig. 3, we show the time evolution of the number 
of condensed atoms, starting with N/Nc = 0.75, found 
by the numerical solution of Eq. (||) with ^ =0.001 and 
7 =0.1, with and without repulsive three-body poten- 
tial. We compare the results of a three-body potential 
with ^3 = 0.016 to the case considered in Ref. |Q, with 
(73 =0. In both, Nc is the critical number for 173 =0. The 
first striking feature with repulsive three-body force is the 
smoothness of the compression mode in comparison with 
the results of 53 = 0. This is a result of the explosive force 
from the repulsion, which oppose to the sudden density 
increase and damps the loss of atoms due to three-body 
recombination effects. Even for lower than 0.016, and 
much closer to (73 = 0, the collapses can no longer "burn" 
the same number of atoms as in the case of 53 = 0. By 
extending our calculation presented in Fig. 3 for all cases 
with > 0.01 and for times beyond ujt = 50, we have 
checked that the number of atoms will increase without 
limit while the condensate is oscillating with frequency 
about 2uj. In particular, the present approach indicates 



that the experimental recent observation of the maximum 
number of ^Li atoms is compatible with 173 much smaller 
than 0.01. The mean square radius for = 0, after 
each strong collapse (when N > Nc) begins to oscillate 
at an increased average radius. The collapse "burns" the 
atoms in the states with higher densities and explain the 
sudden increase of the square radius after each compres- 
sion, remaining the atoms in dilute states. The inclusion 
of the repulsive three-body force, still maintains the os- 
cillatory mode, but the compression is not as dramatic 
as in the former case, and consequently atoms in higher 
density states are not so efficiently burned. The increase 
of the mean square radius (averaged with time) is smaller 
than the one found with only attractive two-body force. 
This is a remarkable feature of the stabilizing effect of the 
repulsive three-body force allowing the presence of states 
with higher densities, as we found in the stationary study. 




FIG. 3. Number of condensed atoms and the corresponding 
mean square radius (r^) (in units of h/{2mu)) as a function 
of time, for ^ =0.001, 7 =0.1. Nc is the maximum number of 
atoms for g-j, = Q. 

Finally, we have to consider that, in the situation when 
the 3-body repulsion dominates over the 2-body attrac- 
tion, the condensate can be in a denser phase where it 
is expected to be strongly unstable due to recombina- 
tion losses. The decay time of the condensate in a denser 
phase is expected to be much smaller than the decay time 
of the condensate in the less dense phase. However, we 
should observe that the dynamics of the condensate is 
modulated by an oscillatory mode with a frequency of 
the order of 2a;, which was already identified by |^8| to 
be ~ even when 53 = 0. In case of g^ > 0, such os- 
cillatory mode dominates the time evolution of the con- 
densate. As the oscillations allow changes in the density, 
the condensate does not "burn" as fast as expected. 

In order to study the condensate decay, we consider 
the original NLSE with the dissipative term and allow 
different possibilities for the three-body interaction 53. 
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We use ^ = 0.001 (the same value used in Ref. @). 
In Fig. 4 we show the result of this study for = 
and Qs = 0.016. The initial number of atoms N can be 
obtained from n, given in the figure. For 53 — 0, we 
took n = 1.625, which is close to the critical limit. For 
(73 — 0.016, we consider three cases: two of them starting 
with the same number of atoms, n — 1.756, but in dif- 
ferent phases (the corresponding chemical potentials are 
P = —1.2, in a denser phase, and (3 — 0.3); and another 
in an even denser phase, with n = 1.965 and /3 = —2.3 
(see also Fig. 1). Based on the results obtained in these 
four different cases, we can estimate that the mean-life 
for the condensate, which is initially in a denser phase, 
is not as small as expected when comparing with gs — 0. 
We observe in this case the relevant role of the oscilla- 
tory mode, related to the frequency of the trap potential, 
which dominates the dynamics of the condensate when 
.93 > 0. 
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FIG. 4. Condensate decay. Number of condensed atoms 
and mean square radius as function of ut, for 7 = and 
^ = 0.001. The chemical potential (3 of the initial state and 
the strengths gs are given in the plot. 

To summarize, our present results can be relevant 
to determine a possible clear signature of the presence 
of a repulsive three-body interaction in Bose condensed 
states. It points out to a new type of phase transition 
between two Bose fluids. Because of the condensation of 
the atoms in a single wave-function this transition may 
present very peculiar fluctuations and correlation prop- 
erties. As a consequence, it may fall into a different uni- 
versality class than the standard liquid-gas phase tran- 
sition, which are strongly affected by many-body corre- 
lations. The characterization of the two-phases through 
their energies, chemical potentials, central densities and 
radius were also given for several values of the three-body 
parameter (73. We develop a scenario of collapse which 
includes both three-body recombination and three-body 
repulsive interaction. From the time-dependent analysis, 
we show that the decay time of the condensate which 



begins in a denser phase is long enough to allow ob- 
servation. However, the observed strongly oscillating 
states are quite different from the analysed stationary 
states. In accordance with the observed strong oscilla- 
tions of the mean-squared-radius, the condensate density 
also strongly oscillates and the observed states cannot be 
characterized as "dense" or "dilute" , justifying the long 
decay time. Nevertheless, through the amplitude of the 
oscillations one can distinguis if the system starts in a 
denser phase. 
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